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Where are the hot ion lines in classical T Tauri stars formed ? 
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ABSTRACT 

Context. Classical T Tauri stars (hereafter CTTS) show a plethora of in- and outflow signatures in a variety of wavelength 
bands. 

Aims. In order to constrain gas velocities and temperatures, we analyse the emission in the hot ion lines. 
Methods. We use all available archival FUSE spectra of CTTS to measure the widths, fluxes and shifts of the detected 
hot ion lines and complement these data with HST/GHRS and HST/STIS data. We present theoretical estimates of 
the temperatures reached in possible emission models such as jets, winds, disks and accretion funnels and look for 
correlations with X-ray lines and absorption properties. 

Results. We find line shifts in the range from -170 km s^^ to -1-100 km s~^. Most linewidths exceed the stellar rotational 
broadening. Those CTTS with blue-shifted lines also show excess absorption in X-rays. CTTS can be distinguished 
from main sequence (hereafter MS) stars by their large ratio of the O vil to O VI luminosities. 

Conclusions. No single emission mechanism can be found for all objects. The properties of those stars with blue-shifted 
lines are compatible with an origin in a shock-heated dust-depleted outflow. 
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1. Introduction 

T Tauri stars (TTS) are young (< 10 Myr), low mass 
(M* < 3Mq), pre-main sequence stars exhibiting strong 
Ha emission. The subclass of Classical TTS (CTTS) are 
accreting material from a surrounding disk. This disk is 
truncated close to the corotation radius by interaction with 
the magnetic field. At the inner rim, the material is ionised 
by the stellar radiation and loaded o nto the field lines 
(jUchida fc Shibatalll984 lKoenigl|[T99ll ) and accelerated to 
nearly free-fall velocity along an accretion funnel before hit- 
ting the stellar surface where an accretion shock forms. 
Angular momentum may be transferred from the star to 
the disk by magnetic torque (|Shu et al.lll994l ). at the same 
time this process is expected to drive out flows possibly in 
the fo rm of a disk wind or bipolar jets. IRomanova et al.l 
(|2004j ) performed full magneto-hydrodynamic (MHD) cal- 
culations of the accretion geometry, while the accretion 
shock and its pos t-shock cooling zone were simulated by 
different authors (ICalvet fc Gullbrind 119981: iLamzinl [l99l 
iGiinther et "alll2007l b 

Observational support for the scenario described above 
is found in data from different wavelength bands. The 
shocked material is heated up to temperatures in the 
MK range and should produce copious amounts of X- 
ray emission. Density-sensiti ve X-ray line ratios strongly 
suggest a post-shock origin (|Giinther et al.l l2007l ) of the 
cool parts of the plasma observed at X-ray wavelengths, 
and the cooling matter veils absorption lin es in the op- 
tical and infrared (jCalvet fc Gullbriii3ll998t ). The geome- 
try of the pre-shock accretion funnel can b e traced by Ha 
line profile analysis (|Muzerolle et al.l |2000D . Undoubtedly 



CTT S also show outflows (a review is given by I Bally et al.l 
|2007[ ). however, the precise origin and the physical driv- 
ing mechanism(s) are uncertain. Theoretical models pro- 
pose a variety of coUimated stellar outflows and disk winds 
to remove angular momentum from the system (IShu et al. 
1994 iMatt et al.l l2002t iGomez de Castro fc Ferro-Fontai] 



20051; IPudritz et all |2007|) . The outflows can take the 
form of winds observed from radio to the UV (e.g. 
Alencar fc Basril 20001: Beristain et al.l 20011 Lanizin et al.l 
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2004 iDupree et al.l l2005al : lEdwards et al.r[20'06[ ). or coL 
limated, often bipo lar jets, again o b served at di f ferent 
wavelength band s (jRodriguezj 119951: ICoffev et all 12004 
iGiidel et al.ll2005l ). The dynamics of the surrounding gas 
have been successfully probed by UV spectroscopy with 
HST/GHRS, HST/STIS and FUSE, especially in the H2 
lines, which show out flows an d fluorescence on the hot disk 
surfa ce (jArdila et all l2002allbl : iHerczeg et all [20 02. 2001 
I2OO6D . In the specific case of TW Hva. lDupree et al. [2005af ) 
claim the existence of a warm/hot wind. Their conclu- 
sion assumes that the CHI and OVI lines have an in- 
trinsic Gaussian shape, with centroid matching the stel- 
lar rest-frame. The observed asymmetry of the CHI and 
OVI line-profiles can then be explained by the presence of 
spherically-symmetric and smoothly-accelerated hot wind 
absorbing the radiation in the blue wing of the se particular 
lines. However, Johns-KruU fc Hcrczc^ (|2007l ) argue that 
this model is incompatible with HST/STIS observations, 
especially for the Civ 1550 A doublet. The existence of a 
hot wind in TW Hya, therefore, remains an open issue. 

In this paper we analyse the shifts and widths of hot-ion 
spectral lines, particularly C ill 977 A and the O vi doublet 
at 1032 A and 1038 A, observed with FUSE. This com- 
plements earlier studies of all pre-rn ain seque nce stars ob- 
served with lUE (jValenti et al.ll2000f ). GHRS (jArdila et al.l 
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2002allH) and of the H2 emission in CTTS ()Herczeg et alj 
20061 ). In Sect. m we briefly summarise the properties of the 
sample stars, and in Sect. [3] we describe the observations 
and data that we use in our present analysis. In Sect. [4] the 
line profile analysis is presented, several possible formation 
regions for the observed properties are discussed in Sect. [3 
and in Sect. [6] we summarise our results. 

2. Stellar properties 

In this section we present the most important properties of 
our sample stars in Table [1] and provide a short summary 
of their main characteristics. 



2.4. V4046Sgr 

ISallmen et al.l (|2000D present the FUSE spectrum of 
V4046 Sgr. V4 046 Sgr is a close b inary system with separa- 
tion of « 9i?Q (|Quast et al.ll2000l ). Both of the components 
are almost identical CTTS, which are accreting matter from 
a circumbinary disk. Chandra X-ray observations show line 
ratios indicative of high density plasma, which cannot be 
explained by coro nal activity, but must be attributed to hot 
accretion shocks (jCiinther et al.|[2006l ). 



2.5. TWA 5 



2.1. RU Lupi 

RU Lupi is a strong accretor, heavi ly-veiled and strongly 
varia ble at different wavelengths (jStempels fc Piskunovl 
|2002() . Its period is uncertain. A wind can be studied us- 
ing cool, metallic lines with a low first-ionisation poten- 
tial, and in II2 as demonstrated by 'H erczeg et al.l (|2005( ) . 
The line profiles of hot ion lin e s in t he UV exhibit com- 
plex structures. iTakami et al.l (|200lh show by spectro- 
astrometry that the Ha emission is significantly extended 
and compatible with a bipolar outflow. In X-rays, RU Lupi 
is observed to be a typic al member of the CTTS class 
(jRobrade fc Schmittll2007f ): speciflcally, it shows a low f/i- 
ratio in the O VII triplet as well as a soft X-ray excess. 



2.2. T Tau 

T Tau is the defining member of the class of TTS. It is a 
hierarchical triple system, where the northern component 
T Tau N is optically visible, while the southern component 
T Ta u S, contributin g mainly in the infrared, is itself a bi- 
nary (jKoreskoll2000D . The T Tau system shows a very sig- 
nificant soft X-ray excess withou t the high densities seen 
in TW Hya (|Gudel et al.l l2007t ). Nevertheless, accretion 
and winds are detectable usin g data for other wavele ngths 
(Ivan Langeve de et allll994bl |al: iDupree et al.ll2005ah . and 
iHerczeg et al.l (|2006[ ). for example, attribute observed fiu- 
orescent H2 emission to outfl ows. The flrs t analy sis of this 
FUSE data, completed by Wil kinson et al.l (|2002l ). provided 
the flrst discovery of extra-solar, H2 Werner-band emission. 



ISallmen et al.l (|2000t ) present the FUSE spectrum of 
TWA 5. TWA 5 is a hierarchical multiple system, the 
main component being a Ml. 5 dwarf. On the one hand, 
TWA 5 has no si gnificant infrared excess indic ating the 
absence of a disk ( Metchev et al.l |2004 [Weinberger et all 



2004: Uchida et al 



20041 ): also, X-ray observations by 



Argiroffi et al.l (|2005l ) suggest that the emission originates 



in a low-density region, not in an accretion shock. On the 
other hand, TWA 5 has an Ha equivalent width of 13.5 A 
and signatures of outflows (jMohantv et al.|[2003( ) , therefore 
we regard it as a transition object. 



2.6. GM Aur 

GM Aur is one of the most popular syst ems to observe stel- 
lar disk constitution and c hemistry (e.g. iBergin et aLll2004l : 
Hueso fc Guil lot 2005; S alvk et al.ll2007D . Numerous obser- 
vations in the radio use the disk of GM Aur to analyse the 
distribution and evolution of dust grains in the disk (e.g. 
iRodmann eralll2006t ISchegerer et aI]|2006D . Based on an 
optically-thi n inner disk surrounded by an optically-thick 
outer region, iNaiita et all (200 7) classify GM Aur as a tran- 
sition object. X-ray observations for GM Aur are discussed 
in Sect. 



2.7. TW Hya 



2.3. DF Tau 

DF Tau is a binary system c onsisting of t wo M stars, re- 
solved for the flrst time by I Chen et al.l (|l990l ). The or- 
bital parame ters are still uncertain after several revisions 
([Schaefer et a l. ,2006) and the distance of the system is still 
a matter of debate. The mo s t reas onable estimate is pro- 
vided bv iBertout fc Genoval (1200^1 . who showed DF Tau 
to be a member of the Taurus-Aurigae moving group 
(w 140 pc) . in spite of its low HIPPARCOS parallax of 
39 ± 7 pc. iLamzin et all (l2001h present HST/STIS and 
lUE spectra showing accretion and wind signatures and 
attribute DF Tau's erratic photometric variability to un- 
steady accretion. DF Tau i s one of few Doppler -mapped 
CTTS (|Unruh et al.lll998h . IHerczeg et aH (|200 6') demon- 
strated that fluorescent H2 emission found in FUSE obser- 
vations is consistent with models of a warm, disk surface. 
X-ray observations for DF Tau are discussed in Sect. 13.31 



Because of its proximity and the absence of a surround- 
ing dark molecular cloud, TW Hya is a key system for 
the study of CTTS. Although relatively old by comparison 
(« 1 Myr), it is still act ively accreting from a surrounding 
disk. Eisner et al. (l2006t) interferometricallv resolve its disk 
and Hcrczcg et al.l (|2004| ) constrain the inner disk temper- 
ature to about 2500 K using H2 fluorescence in HST/STIS 
and FUSE data. Observations with STIS and FUSE show, 
in addition to H2 , atomic emission lines with notably asym- 
metric shapes. The analytical models expect a co ld wind 
to be driven by magnetic flelds from the disk (iShu et al.l 
Il994[ ). and this can be tr aced observationally from sin- 
gle or double ionised ions (jLamzin et al.l [20041 1. TW Hya 
has been extensively observed in X-rays, which revealed a 
comparatively cool and dense emissio n region i nterpreted 
as signature of an accre tion hot spot (jKastner e t al. 2003; 
IStelzer fc SchmittI I2004D. Our own theoretical modelling 
supports this view (iGiinther et~ani2007l ). 
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Star 


Spectral type 


d 


'^rad 


V sm I 




Nh 


Av 






[pc] 


n — 1 1 

[km s ] 


[km s ] 


[°] 


[10 cm ] 


[mag] 


RU Lup 


G5V 


140 


-1.9 ±0.02 (17) 


9.1 ±0.9 (17) 


24 (19) 


1.8 (15) 


0.07 (9) 


T Tau 


G5V 


140 


18 (9) 


20 ± 5 (9) 


15 (9) 


5 (5) 


0.3 (9) 


DF Tau 


MO-2 


140 


21.6 ±9.2 (4) 


15.8 ± 1.0 (7) 


> 60 (21) 


3-8 


0.5 (9) 


V4046 Sgr 


K5 


83 


-6.94 (12) 


14.9 ±0.9 (18) 


35 (18) 


<2 (6) 


0.0 (18) 


TWA 5 


M3V 


57 


14 (20) 


uncertain (20) 


n.a. 


0.3 (2) 


n.a. 


GM Aur 


K5V 


140 


24 (10) 


13 (3) 


56 ±2 (16) 


0.2-0.4 


0.31 (11) 


TW Hya 


K8V 


57 


12.5 ±0.5 (1) 


5±2 (1) 


7± 1 (13) 


0.35 (14) 


0.0 (8) 



Refer ences fP lAlencar fc Batalliall2002l: (2) lAjg iroffi ct al. 20 051: (3)lBas ri & Bataiha 1990; ( 4)lBouvier et al.lll986l : ( 5)lGiidel et al 



20071 : (6) new f its to the data from 'Giinther et al.,.2006: ('7') .Hartmann fc St auffcr 1989; f8VH erczeg et al.ll200 jrf9H Herc zeg et al 



200d: (10) llovl 119491: (11 ) ||luzcrollc ct al. 1991(12) iQuast et al.ll2000l: (13 ) Qi ct al. 2004; ( 14) lllobrade fc Schm itt 20 ol; (15) 



Robrade fc Schmittl2007l: (le i iSimon et al.ii200Cr ('17}IStemDels fc Piskunwr2002i : (;i8i iStempels fc Gahml2004l : Cigl iStempeis et all 
2001 : (20) [Torres et al.ll2003 : (21) lUnruh et al.lll998l 



3. Observations 

3.1. FUSE data 

In order to perform a systematic study of the FUV proper- 
ties of CTTS, we retrieved all available FUSE observations 
of TTS from the archive. For a comparison with "normal" 
stars, we selected two main sequence stars, the rapidly- 
rotating K star AB Dor and the solar analog a Cen A. 
TableOlists the observation logs; all targets except a Cen A 
are observed using the LWRS aperture with 30" x 30" 
field of view We reduced the data with CalFUSE v3.2.0 
(|Dixon et al.ll2007l ). performing extractions of the complete 
data set and night-only time intervals to check for airglow 
contamination. With the exception of GM Aur, the spectral 
regions of interest are not affected by airglow. To coadd the 
individual exposures, we used the CalFUSE get_shif t task 
and applied two different methods: both methods involved 
cross-correlation, the first using a stellar emission line (O vi 
1032 A on the laLiF detector and C ill 977 A on the 2bSiC 
detector), and the second method using a strong airglow 
line (H Ly/3 and H Ly7). Airglow lines completely cover 
the aperture and cannot therefore be used to judge cen- 
tering of a target. We studied all stellar spectra with OVI 
and CIII sufficiently strong for cross-correlation, and found 
no significant differences between the line profiles derived 
using both methods. Using airglow lines does, however, at 
times enable H2 emission to be detected in the IbLiF chan- 
nel. This is in contrast to the less accurate cross-correlation 
achieved using low S/N stellar lines. 

It seems that in our observations the targets did not 
drift through the apertures. Therefore we shift all airglow 
H Ly/3 lines to match the laboratory wavelength. As the 
targets may be positioned differe ntly iri aperture we per- 
form the following check: Herczeg 'et al.l ((^02) discovered 
molecular hydrogen at rest compared to the star, so we 
identify H2 lines in the IbLiF channel in TW Hya (the 
shift between ch annels a and b on t he same chip is smaller 
than 8 km s"! ijPixon et al.ll2007D ): the shift of our H2 
lines should now match the radial velocity of TW Hya. 
This fixes our wavelength scale and from here on we use 
TW Hya as our reference spectrum. To avoid the need for 
a spectral identification, we search the IbLiF channel of the 
other observations for strong features, which correlate with 
TW Hya. Only lines with the same strengths in full and 
night-only extractions are considered here. We also check 
that no airglow feature is expected at that position accord- 



ing to iFeldman et all (12001). The line wi dths we find are 
as small as expected (jHerczeg et al.|[2006[) . lending further 
credibility to the coaddition of single exposures using an 
airglow line. In addition, we cross-correlate larger regions 
of the spectra without strong airglow features and find that 
the required shifts are in agreement with the analysis of 
single lines. For RU Lup, T Tau and DF Tau we correct 
for the velocity of th e H2 outfiow relative to the star from 
iHerczeg et al.l (j2006l) . Taking into account the stellar radial 
velocities, the lines shifts confirm the shifts adopted above 
from the H Ly(3 lines. For V4046 Sgr, we do not have a 
HST/STIS measurement of H2 lines, and for TWA 5 and 
GM Aur no suitable features in the IbLiF channel were 
found. However, given the success of the wavelength cali- 
bration for other stars, we are quite confident that it is re- 
liable also in these cases. We estimate the total wavelength 
error in the la/bLiF channels to be 20 km s~^. We then 
correlate the 2bSiC channel with the laLiF channel using 
the stellar Ovi 1032 A line. Unfortunately the signal-to- 
noise ratio (hereafter SNR) in 2bSiC is much lower than in 
laLiF and for T Tau, DF Tau and GM Aur we have to use 
the H Ly/3 airglow line. Because the FUSE channels are 
not perfectly aligned, the position of the target in the vari- 
ous channels might be different; the wavelength calibration 
below 1000 A for these three stars is therefore uncertain. 

For data of low SNR CalFUSE sometimes overcorrects 
the background, leading to negative flux intensity values; 
we therefore always separately fit the background in line- 
free regions, close to the line of interest. 

3.2. HST data 

Four of our sample stars (RU Lup, T Tau, DF Tau 
and TW Hya) have been observed with gratings onboard 
of the HST. RU Lup, T Tau and DF Tau were ob- 
served with the GHRS; these observations are p resented 
by lArdila et all (|2002af) and iLamzin et all (|2001l ) and we 
adopt the linewidths and shifts measured by these authors. 
The STIS observations of RU Lup were carefully anal- 
vsed bv iHerczeg et all (120 Q5D and those of TW Hya b y 



IHerczeg et al.l (|2002[) and iJohns-Krull fc Herczed (|2007D . 



We use the STI S spec tra of T Tau and DF Tau from 
CoolCAT (lA^[2005l) . which is a database of reduced 
and calibrated UV spectra taken with STIS for cool stars. 
These spectra all belong to the HST program 8157 and 
were taken in the year 2000. The observations of TW Hya 
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Table 2. FUSE observations 



Target 


Date 


Exposure time 


Data ID 






[ksec] 




RU Lup 


2001-08-28 


24 


A1090202 


T Tau 


2001-01-15 


21 


P1630101 


DF Tau 


2000-09-13 


26 


A1090101 


V4046 Sgr 


2000-05-18 


16 


P1920202 


TWA 5 


2000-05-15 


15 


P1920101 


GM Aur 


2004-03-03 


32 


D0810101 


GM Aur 


2004-10-23 


34 


D0810102 


1 W -tiya 


2000-06-03 


2 


rl860101 


TW Hya 


2003-02-20 


15 


G0670101 


TW Hya 


2003-02-21 


15 


G0670102 


AB Dor 


1999-10-20 


22 


X0250201 


AB Dor 


1999-12-14 


24 


X0250203 


AB Dor 


2003-12-26 


102 


D1260101 


a Gen A^ 


2001-06-25 


15 


P1042601 


a Gen A^ 


2006-05-05 


13 


G0810102 



Notes: (1) MDRS aperture 



are not contained in CoolCAT, because they were taken 
in a non-standard setup. For a consistency check, we take 
T Tau, and find that the spectra pr esented there agree w ith 
the previously pubhshed results in iHerczeg et al.l (|2005l ) . 

3.3. X-ray data 

Most of our sample stars have been observed at X-ray wave- 
lengths with Chandra, XMM-Newton or both, and we ob- 
tained these data from the literature. For DF Tau and 
GM Aur, there are no X-ray spectra available in the lit- 
erature, however, both DF Tau and GM Aur have been 
observed with ROSAT/PSPC pointings. We therefore re- 
trieved the ROSAT X-ray data from the archive in order 
to enlarge our sample for comparisons of optical and X-ray 
absorption. In Table [31 we provide information on the data 
used. The ROSAT data were reduced using EXSAS version 
O30ct. Source photons were extracted from a circle with a 
radius of 2.5 arcmin to avoid contamination from nearby 
sources. The background was taken from a region close 
to the target without detected sources with a radius three 
times larger. We measured fluxes of 0.013 ± 0.002 cts s~^ 
for DF Tau and 0.021 ± 0.003 cts s'^ for GM Aur. Using 
MIDAS/EXSAS, we fitted Raymond-Smith and Mewe- 
Kaastra models assuming cold absorption. The ROSAT 
data are of low SNR. We therefore tried to use different 
data binnings and compared our final results. We provide 
a range of possible absorption column-densities in Table [1] 
because, as often found, observationally a large amount of 
cool plasma and small interstellar absorption lead to very 
similar signatures. To test the reliability of our procedure, 
we fit the ROSAT data for AB Aur, which is located in the 
same field as GM Aur. We measure an absorption column- 
denstiy Nh of 2 — 3 x 10^° cm~^, which is approximately 
half of the m easurement der i ved u sing XMM-Newton data 
presented bv lTelleschi et all (|2007[ ). We therefore conclude 
that, taking into account systematic errors, the ROSAT- 
derived absorption columns are consistent with those de- 
rived using XMM-Newton data. 



Table 3. ROSAT observations 



Target 


Date 


Exposure time 


Obs ID 






[ksec] 




DF Tau 


1993-08-93 


10 


WG201533P.N1 


GM Aur 


1992-09-22 


1.3 


WG201278P.N1 


GM Aur 


1993-03-04 


4 


WG201278P-1.N1 



4. Results 

4.1. FUV lines 

In Fig. [H we present the C ill line profiles at 977 A and 
both components of the O vi doublet at 1032 A and 1038 A 
for our targeted stars. Each spectral line is fitted with a 
Gaussian profile and shown normalised with respect to the 
fitted peak value. The best-fit parameters for line flux, shift 
of the Gaussian-line centre, and its full width at half max- 
imum (FWHM) , ordered by the observed shift in the O vi 
1032 A fine, are given in Table [H where the shift is cor- 
rected for the radial velocity of the star. The statistical 
errors on the shift and the FWHM given in the table are 
typically comparable to our systematic calibration uncer- 
tainties. The systematic errors on the fluxes are of the order 
10-20%. The Ovi doublet line at 1032 A is detected in all 
sample stars and the C ill line in all except GM Aur. In all 
cases, the line shifts measured in the C III and O vi are con- 
sistent, arguing for related emission regions; only for T Tau 
is there a statistically significant difference, probably due to 
uncertainties in the wavelength calibration of 2SiC detector 
data, as discussed in Section [3T] 

In the case of RU Lupi, an absorption component is 
clearly visible in all three lines, which is also fitted with a 
Gaussian absorption line. We note that this absorption is 
not self-absorption. [Herczeg et al.l (j2005( ) identify the ab- 
sorption in the Ovi 1032 A as H2 6-0 P(3) shifted by 
—23 ± 5 km s~^ relative to the H2 rest wavelength, and 
the absorption i n O yi 10 38 A as Cii. Using the package 
H2oolts (, McCandhssl2003t ) . we identify the absorption com- 
ponent in C III to be a superposition of the equally strong 
H2 2-0 P(5) and 11-0 R(3) shifted by « -10 ± 15 km s"! 
again relative to the H2 rest wavelength. Within the Ovi 
1038 A line of TW Hya, we fit the C 11 blend on the blue side 
with an additional Gaussian component. Unfortunately the 
data quality does not enable further blends and features to 
be identified. In the FUSE bandpass, numerous fiuorescent 
Lyman and Werner ban d H2 lines are prese nt, which are 
usually relatively sharp ( Herczeg et al.|[2006 ). The lines in 
Fig. [1] are well-described by one-peaked distributions, so 
the contribution of any individual blending line has to be 
small compared to the broad ion lines. 

Atomic physics predicts that the 1032 A component of 
the Ovi doublet should be twice as strong as the 1038 A 
component. We can confirm this prediction for the data 
of V4046 Sgr and TWA 5. In Table H significant devia- 
tions from the expected line ratio can be seen in RU Lup 
and TW Hya, where additional components complicate the 
fit. However, in the cases of T Tau, DF Tau and GM Aur 
the Ovi 1038 A line is almost completely absent, because 
of absorption by larger molecular hydrogen columns (see 
TableU]) in the lines H2 5-0 R(l) at 1037.1 A and 5-0 P(l) at 
1038.2 A. RU Lup and T Tau, both CTTS with blueshifted 
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Table 4. FUSE best fit parameters, fluxes in 10^^^ erg cm"^ s~^, widtli and sliift in km s^^, errors are statistical only 







Cm 977 A 






Ovi 1032 A 






Ovi 1038 A 




Cm 1175 A 


Star 


Flux 


Shift 


FWHM 


Flux 


Shift 


FWHM 


Flux 


Shift 


FWHM 


Flux 


RTT T iin 


1 U 31 UU 


1 9n + 1 n 




Z4; ZIZ U 


1 7n + 9n 

LIU ZL ZU 


Qnn + fin 


1 o zn o 


1 nn + 90 

— lUU ZIZ zu 


1 / U ZIZ ou 


9n + 9 

ZU ZIZ z 


- Absorption 


-60 ± 40 


-130 ± 10 


100 ± 20 


-4±0 


-200 ± 6 


50 ±20 


-10 ±2 


-100 ± 10 


120 ± 30 




T Tau 


10 ±4 


1 ± 13 


90 ±30 


8±2 


-50 ± 10 


130 ± 30 


n.a. 


n.a. 


n.a. 


18 ±2 


DF Tau 


17 ±6 


-6± 12 


120 ± 30 


11±3 


-22 ± 16 


200 ± 40 


n.a. 


n.a. 


n.a. 


44 ±2 


V4046 Sgr 


60 ±20 


50 ±20 


310 ± 70 


37 ±6 


19 ±8 


160 ± 20 


16 ±4 


1 ± 12 


150 ± 30 


50 ±3 


TWA 5 


30 ±10 


60 ±20 


160 ± 40 


30 ±5 


30 ± 7 


130 ± 20 


14 ±3 


19 ±9 


110 ±20 


25 ± 2 


GM Aur 


n.a. 


n.a. 


n.a. 


5± 1 


60 ±20 


210 ± 60 


n.a. 


n.a. 


n.a. 


24 ±2 


TW Hya 


200 ± 25 


100 ±6 


190 ± 20 


220 ± 20 


74 ±5 


220 ± 15 


100 ± 15 


35± 15 


290 ± 40 


310 ±4 


- Blend 














9±9 


-140 ± 14 


80 ±40 





VI 1032 
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RU Lupi 
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-10 L 
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lombdo [A] lambda [A] 
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lombdo [A] lambda [A] 



TW Hya 




-400 -200 200 400 -400 -200 200 400 -400 -200 200 400 
V [km/s] V [km/s] v [km/s] 



1 173 1 174 1 175 1 176 1 177 1 1 73 1 179 
lambda [A] 



Fig. 1. Hot ion lines observed with FUSE in CTTS and 
best fit Gaussian profiles ordered by the shift in the Ovi 
1032 A line. All hne profiles are normalised to tlieir the 
Gaussian profile peak value and rebinned to instrumental 
resolution. 



emission, are found to have the largest reddening and ab- 
sorption column density. 

In Fig. [51 we present the FUSE observations of the C ill 
multiplet at wavelength 1175 A. Unfortunately the SNR of 
the data is too low to be able to resolve the H2 lines, and 
to fit the individual line components. This analysis would 



Fig. 2. The Cm multiplet around 1175 A rebinned to in- 
strumental resolution. 



have enabled density diagnostics to be measured. For RU 
Lup and T Tau in fact, the multiplet is barely detectable. 
Therefore a detailed Gaussian function fit yields no useful 
information, and in Table [4] we present only the total flux 
and its statistical error. The dereddened ratio of the fluxes 
in the C ill 977 A line and the 1175 A multiplet is, in princi- 
ple, density-sensitive in the r ange 8 < log n < 11 according 
to th e CHIANTI database ()Dere et al.l 119981: iLandi et al.l 
|2006[ ). However, apart from RU Lup the observed values are 
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-400 -200 200 400 -400 -200 200 400 -400 -200 200 400 
. [km/s] . [km/s] . [k,T,/s] 

Fig. 3. Hot ion lines in the STIS spectra. The relative flux 
of the red doublet member is multiplied with a factor of 2 
(red/grey hues). 



far below the predicted ratios. iRedfield et al.1 (|2002[ ) have 
analysed this line ratio in FUS E observations of late-type 
dwarf stars, and iDupree et al.l ([2005b) for luminous cool 
stars. Both authors conclude that the Cm 977 A line is 
optically-thick in the transition region and that in addition 
there is strong interstellar line absorption. 

The STIS data from CoolCAT is shown close to the hot 
ionic doublets Nv 1240 A, Siiv 1400 A and Cv 1550 A in 
Fig. [21 In each case, the expected intensity ratio is two 
to one. To facilitate the comparison of the line profiles we 
show the flux of the red doublet member increased by a 
factor of two. In the absence of blending and interstellar 
absorption, both lines should show a perfect match. The 
red member of the Nv doublet at 1242.8 A is totally ab- 
sorbed by an interstellar Ni line, so we do not attempt to 
fit the line profile. For the other lines, the results of a fit 
with a single gaussian component are shown in Table [51 In 
DF Tau there is no overlap between the two lines shown in 
the middle panel of Fig. [31 Therefore they cannot be dou- 
blet members, but have to be other lines, likely H2. These 
lines also distort the red member of this doublet in the 
other two stars, leading to a shifted centre in the gaussian 
fit. The large errors result from generally non-Gaussian line 
shapes. The Table also includes the GHRS data and the ob- 
servations of TW Hya, resembling the FUSE data. These 
profiles are nearly triangular, rising sharply close to the rest 
wavelength and ex1 :ending to about 400 km s~^ on the red 
side (see Fig. 7 in iHerczeg et al.l [2002D . Of similar shape, 
although extending only to 200 km s~ ^ on the red side, is 
the Civ 1550 A doublet in DF Tau (jLamzin et al.ll200l[ ) 
in the Gfffl^' observation. This da ta is fitted with multiple 
Gaussians bv lArdila et al.l (|2002al ). We give here the values 
for the dominant component only. 

Despite the uncertainties in the profiles of individual 
lines, it is obvious from Tables [H and [5l and Fig. [1] and [H 
that FUV lines in CTTS cover a range of wavelength shifts 
from blue-shifted to red-shifted velocities of the order of 
-200 km out to 100 km s~^, and could therefore poten- 
tially originate in in- or outflows or both. RU Lup clearly 
shows blue-shifted emission, the same is likely the case for 
T Tau, because the O vi shift can be more accurately de- 



termined than the Cm shift. DF Tau and V4046 Sgr are 
consistent with centred lines and TWA 5, GM Aur and 
TW Hya with red-shifted lines. 

4.2. Excess absorption 

In Fig. m we show the absorbing column densities as fltted 
to the X-ray observations and optical Ay values. With a 
standard gas-to-dust ratio, they should b e related through 
the formula Nh = Ay ■ 2 x lO ^^cm"^ (jSavage fc Mathid 
Il979i . but see lVuong et al.l (|2003D for a compilation of other 
conversion factors in the literature, all roughly consistent 
with this value). The determination of Ay values is diffi- 
cult for CTTS, where not only the spectral type, but also 
the veiling needs to be known to compute the intrinsic 
colours. This renders Ay values in the literature notori- 
ously uncertain, especially if they are based only on pho- 
tometric and not on spectroscopic information. We have 
tried to compile the best available estimates in Table [H 
For AA Tau and possibly other CTTS, the circumstellar 
absorption and the veiling are time-d ependent. In Fig. [H 
we sho w m ean values for AA Ta u from lSchmitt &: Robradd 
(|2007D and iBouvier et all (|1999D . The Ay for T Tau is of- 
ten given one magn itude larger than the value we use, but 
[Walter et al.l (|2003l ) convincingly show that this is an over- 
estim ate. In the case of DF Tau, most values are between 
0.21 ([Kenvon fc HartmannI I1995D and 0.5 5 (IStrom et al 



1989), with the exception of 1.9 ± 0.6 in ICohen fc Kuhi 



(1979): for RU Lup, all mo dern mea s urements agree on the 
low Ay value dL amzin c t'al] [19961 : [Stempels fc Piskunovl 
[2002[ : [Herczeg et al.n2005i) . In addition to the TTS stud- 
ied in this paper, we incl uded data for the young stars 
BP Tau and SU A ur from [Robrade fc SchmittI (|2006^ and 
[Errico et al] ([200lh. for BP Tau, AB Aur and HD 163296 
from Telleschi et al. ([20071 see refere nces therein) and for 
MP M us from [Argiroffi et al.[ Hool) and [Mamaiek et al.[ 
([2002h . In the case of AA Tau, reddening and absorb- 
ing column are exceptionally large, because this object is 
seen almost edge-on. AA Tau, DF Tau, T Tau, RU Lup 
and V4046 Sgr all show much larger X-ray absorbing 
columns than expected from the optical extinction, indi- 
cating large amounts of gas with a very small dust content 
compared to the interstellar medium. This effect is signifi- 
cantly larger than the uncertainties in the Ay and Nh val- 
ues. In V4046 Sgr the relative uncertainties are higher com- 
pared to the other stars, because Chandra is not as sensitive 
as XMM-Newton at soft X-ray energies. In RU Lup and, 
possibly, T Tau the hot ion lines exhibit significant blue 
shifts. Therefore, it may well be that the hot UV-emitting 
plasma provides the observed X-ray excess absorption. 

All Nh values provided are obtained by fitting a cold 
absorber model, implying that the elements contributing 
to the X-ray absorption are not fully ionised. The most 
important elements are O, H, and He; with increasing tem- 
perature and ionisation an even larger column density is 
necessary to provide the observed absorption. 



4.3. Ratio of Ovii to Ovi luminosity 

X-ray data are available for all CTTS sample stars, most of 
which were observed using gratings. The resolution of these 
X-ray lines is insufficient to determine velocities, and we 
therefore analyse only the measured line fluxes. In Fig. [51 we 
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Nv 1239 A 


Si IV 1394 A 


Si IV 1403 A 


C IV 1548 A 


Civ 1551 A 


Star 
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Shift 


FWHM 


Shift 


FWHM 


Shift 


FWHM 


Shift 


FWHM 


Shift 


FWHM 






km 


km s~^ 


km 


km 


km 


km 


km 


km s~^ 


km s^^ 


km s^^ 


RU Lup 


GHRS 






5±5 


270 ± 10 


-31 ±6 


250 ± 10 


-20± 20 


340 ± 40 


-10± 10 


300 ± 30 


RU Lup 


STIS 


-15 ±25 


190 ± 70 


26 ± 20 


200 ± 50 


-20 ± 20 


180 ± 50 


-25 ± 25 


230 ± 80 


15 ±18 


120 ± 50 


T Tau 


GHRS 






-15 ±7 


220 ± 20 


-50 ± 10 


220 ± 30 


-25 ±6 


250 ± 10 
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150 ± 10 
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-45 ± 30 


220 ± 80 
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150 ± 50 
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320 ± 20 
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30 ±40 
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-40 ± 40 


230 ± 100 


-50 ± 40 


180 ± 120 
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redshifted, see text 


42 ±3 


170 ±6 


9±3 


130 ±8 


redshifted, non-g 


aussian, see text 
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Fig. 5. Ratio of dereddened fluxes in the Ovii triplet at 
22 A to the Ovi 1032 A line vs. total luminosity in those 
lines for CTTS (triangles) and other stars (squares). See 
text for data sources. 



Fig. 4. Optical reddening and absorbing column density 
from X-rays. Error bars represent 90% confidence intervals 
on the Nh fit, except for DF Tau and GM Aur where we 
give the (tighter) ranges discussed in Sect. 13.31 The dotted 
line indicates the values predicted using Nh = Ay ■ 2 x 
lO^^cm-^. 



compare the fluxes of the O VI and O Vil lines, where the X- 
ray data is taken from the same refer ences as Nh in Table d 
The observed fluxes are dereddened (jCardelli et al.llT989f ) in 
the UV with a standard dust grain distribution assuming 
Rv = 3.1, and in the X-rays following the col d abso rption 
model by iBalucinska-Church fc McCammonI (|1992|) . For 
comparison, some star s with " n ormal " coronal X-ray emis- 
sion have been added (|Ness et al. 2004 ^, whose O vi emis - 
sion is known dR edfield ejalJlMa; iPupree et al.ll2005tf ). 
Where iNess et al . (2004) provide measurements from both 
XMM-Newton and Chandra, we show two connected sym- 
bols to give an idea of the intrinsic variability of the 
sources. This by far dominates the observational uncer- 
tainties, which are approximately 10%. As is clear from 
Fig. O the CTTS separate well out in the upper right cor- 
ner of the plot, indicatin g an Ovii excess compare d to 
dwarf and giant cool stars iRobrade fc Schmitt) (|2007( ) and 



iGiidel fc Telleschil (|2007h show a very similar analysis us- 
ing the O viii/0 VII ratio again with an O vii excess. These 
studies clearly show that the observed excess radiation of 
CTTS is confined to a relatively narrow temperature range 
about the formation of the He-like O vii triplet at 1-2 MK. 
We note that the CTTS T Tau has an unusually strong 
Ovii excess, with a ratio close to 55. Estimating the H2 
column density to be roughly a quarter of the total neu- 
tral hydrogen absorbing column (|Bloemenlll987[ ). no more 
then about a third of the O vi flux in T Tau could be ab- 
sorbed and the true blueshift cannot exceed 100 km s~^. 
Otherwise emission blueward of the absorbing H2 P(3)(6- 
0) line at 1031.19 A would be observed in Fig.[TJ Absorption 
in the O vi lines cannot therefore explain the extreme O vii 
to O VI ratio. A possible origin for the large emission mea- 
sure at 1-2 MK might be a small-density accretion shock 
with a large filling factor. 



5. Interpretation 

In the following we attempt to constrain the emission re- 
gions and set limits and constraints on possible models for 
the origin of the observed line shapes. 



8 



H.M. Giinther and J.H.M.M. Schmitt: Formation of hot ion lines in CTTS 



5.1. Blue-shifted emission 

The detected emission of hot-ion hnes in the FUSE data 
appears blue-shifted in RU Lupi and, probably, in T Tau. 
Both systems show small values of v sin i compared to zero- 
age main-sequence stars and their inclination is close to 
face-on (Table [Ij. However, TW Hya is a counter-example 
with small inclination and red-shifted emission. 



5.1.1. Shocks in outflows 

We assume that an outflow exists, produced by either a 
disk or stellar wind. Blue-shifted emission could then be 
due to either internal shock fronts or to termination shocks 
of outflows running e.g. into an ISM cloud. For the case 
of RU Lupi, the measured centroid velocity is found to be 
Vohs ~ 200 km s~^ away from the star, and the pe ak for- 
matio n temperature of Ovi is 300 000 K (Mazzot ta et al.l 
I1998D . If we assume th at the O vi gas is hea t ed by a strong 
shock wave, following IZel'Dovich fc Raizeil (|1967f) we can 
infer that the shock front properties from the Rankine- 
Hugoniot jump conditions. The emission of a post-shock 
cooling zone would be observed with velocity Wobs, which 
sets a lower boundary on the shock velocity Ushock to 



4Wobs = V: 



shock 



> 800km s~^ . 



(1) 



By using line shifts we measure only the line of sight 
velocity, but RU Lupi and T Tau are probably viewed 
close to pole-on, thus we expect this value to be close 
to the true velocity. Typ ical jets reach u p to 4 00 km 
(|EislofFel fc Mundtlll998l) . but iRodrieued (|1995D report ve- 
locities of up to 1000 km s~^. The highest bulk velocities ob- 
served in the jet of the CTTS DG Tau are about 500 km s'^ 
(iLavallev-Fouauet et al. 2000). Using the strong shock for- 
mula the post-shock temperature Tpost can be obtained by 



-^post — 



7-1 



(7+1)^ 



2/xmH 2 

1 ^shock 



(2) 



k 



where mn is the mass of an hydrogen atom, /i the dimen- 
sionless mean particle mass, k is Boltzmann's constant and 
7 is the adiabatic index. In the last step, we use 7 = 5/3 for 
an ideal gas and Eqn. [TJ The estimated post-shock temper- 
atures of « 2 X 10^ K for RU Lup and « 2 x 10^ K for T Tau 
are large eno ugh to form O VI an d O vii. This seems plausi- 
ble, because iGiidel et al.l (|2005D present already spatially- 
resolved X-ray emission from a CTTS jet in DG Tau A. 
The shifts of the O vi lines observed in RU Lup and T Tau 
are larger than measured for the C ill lines, agreeing with 
models of post-shock cooling. C ill is formed at lower tem- 
peratures, when the post-shock flow has lost energy, cooled 
and slowed down. 

The observed volume emission measures (in units of 
10^1 cm-3) for RU Lup, T Tau and DF Tau in the Ovi 
1032 A and 1038 A lines are 1.3, 1.4 and 4.7, assum- 
ing the temperat ure of maximum emi ssivity and deredden- 
ing of fluxes by ICardelli et al.l (|1989D : for Cm the corre- 
sponding numbers are 3.9, 1.5 and 8.0. Along the jet of 
the CTTS DG Tau pre-shock densities between 10^ and 
10^ cm~ ^ are observed, with the hig her values closer to 



distance d and we find d « 3 x 10^^ cm for RU Lup with 
V = 800 km s~^ and n = 10'* cm~'^. The diameter of the 
jet is then of the order 10*'^ cm and the column density of 
the cooling zone alone Nh = 3 x 10^" cm~^, adding to the 
absorbing column in the base of the jet, the circumstellar 
material further out and the interstellar column density. 
Therefore the large column densities observed for the stars 
with blue-shifted emission might be partly due to jet mate- 
rial. In T Tau and DF Tau the observed velocities are lower, 
which results in d r:^ 3 x 10^'' cm. This is too small to ex- 
plain the observed emission measures. In the jet emission, 
picture this suggests jets inclined to the line of sight, so the 
velocities and therefore the cooling distances are underesti- 
mated. According to the data in Table [U the inclination of 
T Tau is smaller than the inclination of RU Lup, but it is 
not clear if jets are always orientated perpendicular to the 
accretion disk. DF Tau has the largest inclination, so we 
expect to see the jet at a high inclination. This might ex- 
plain the large emission measure despite the low observed 
shift, which is even consistent with centred emission. 

It could thus well be that the additional absorption of 
the stellar X-ray emission is due to the hot gas from which 
the FUV lines originate. 

5.2. Anomalous line widths 

A comparison of rotational velocities listed in Table [1] 
with the line widths shown in Table [3] and [5] shows that 
the observed lines are superrotationally-broadened despite 
the measuremen t uncertainties. Th is is not uncommon in 
late- type stars (jWood et al.l Il997f ). To illustrate this is- 
sue further, in Fig. [S] we compare the O vi 1032 A line 
measured for TW Hya, the much more rapidly rotat- 
i ng K dwar f AB Dor with vsini = 100 ± 5 km s~* 
(IVilhu et al.1 Il987[ ). and t he slow-rotator a Cen with 
vsini = 2.7 ± 0.7 km s~^ (ISaar fc Ostenlll997D is in fact 



the star (iLavallev-Fouauet et al.1 gOOOf ). We now use the fit- 
formula from lHartigan et al.l ( 19871 ) to estimate the cooling 



comparable to TW Hya. lAke et al.l (|2000r found signif- 
icantly superrotationally-broadened coronal lines already 
in AB Dor and speculate that this might be due to mi- 
croflare heating. The blue wing of TW Hya matches the 
data from AB Dor quite well, but the asymmetric profile 
from TW Hya is significantly different on the red side. At 
any rate, the O vi 1032 A is clearly far broader than the line 
from a Cen. In CTTS, redward line asymmetries are usually 
attributed to emission in the accretion process. We there- 
fore consider if this interpretation is also valid for TW Hya. 



5.3. Magnetospheric infall 
5.3.1. The basic model 

In the standard accretion scenario for CTTS, the disk is 
truncated at a few stellar radii and the temp eratures on the 
disk surface are a few thousand K at most (|Herczeg et al.l 
2004). Consequently, excluding additional energy deposi- 
tion by magnetic star-disk interaction, the accreting mate- 
rial reaches the formation temperature of O VI for the first 
time shortly before it passes through the accretion shock. 
At this stage, the material on the opposite side of the star 
ought to be blocked from view and only the low velocity, 
low temperature zones of the accretion funnel, close to their 
origin on the disk, can be seen blue-shifted if they face the 
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TW Hyo 
AB Dor 
A Cen A 






. [Krn/s] 

Fig. 6. Comparison of the Ovi 1031.91 A line for three 
stars. 



observer. Therefore hot ion lines can be observed only red- 
shifted from the funnel flow. 

Using th e CHIANTI 5.2 database (|Dere et all Il998t 
iLandi et al.l l2006') we have numer ically modelled 1 :he ac cre- 
tion spot as described in det ail by Giinther et al.l (I2007D for 
TW Hya and for V4046 Sgr (|Gunther fc Schmittll2007ir . We 
measure the infall velocities in the accretion funnel to be 
525 km s~^, and the pre-shock densities around 10^^ cm~^ 



-200 200 
velocity shift [l<m/s] 



faO'O 



Fig. 7. O VI line at 1032 A. solid hne: FUSE data, dashed 
line: simulation with purely thermal broadening, dotted 
line: simulation with additional turbulent broadening of 
150 km s-^ 



locities, the post-shock temperatures and hence the emit- 
ted X-ray fluxes are smaller. Such effects are not included in 
our models, but could contribute to the observed FUV flux. 
We are unable to reliably constrain fits of inhomogeneous 
spots due to insufficient diagnostics. 



for TW Hya and 2 x 10" cm for V4046 Sgr. This mat- 5,3^2. Predicted UV line shifts 



ter passes through a shock on the stellar surface and is 
heated up to temperatures of « 2 x 10^ K. X-ray emission 
is produced in the following cooling fiow in the post-shock 
zone. We explicitly calculate the ionisation and recombina- 
tion rates, to be able to assess the ionisation state at each 
depth. Fitting this model to high-resolution X-ray data, we 
measure the abundances of the most important elements 
and determine mass accretion rates and filling factors be- 
tween 0.1%-0.4%. 

Our model explicitly assumes optically thin-radiation 
loss. Along the direction of inf ailing flow, most resonance 
lines are found to be optically-thick up to optical depths 
above 10. The opacity is measured using a simple opacity 
calculation, and its precise value depends on both the infall 
density and elemental abundance. In the He-like triplets of 
Neix and Ovii, observed in X-rays, the ratio of resonance 
line to intercombination plus forbidden line is expected to 
be close to unity from atomic physics. This is matched by 
observations, so we can infer that the resonance lines do not 
suffer much absorption. We note that the detailed accretion 
region geometry is unknown. In particular, the system ge- 
ometry is probably not 1-D and photons may well escape 
through the boundaries of the system rather than along the 
direction of flow. Unfortunately, this leads to considerable 
uncertainty in the computed flux. In our simulation how- 
ever, emission is generated in the post-shock region to the 
same extent as observed for the OVI triplet, but the post- 
shock contribution to the CHI line is small. We caution that 
this conclusion is based on the assumption of very small and 
hot spots. Theoretical si mulations of the infall geometry 
using either dipolar fields (jRomanova eiaD'2004) or more 
realis tic geometries scaled from other stars (Grcgor v et al.l 
l2006t ) suggest inhomogeneous spots with a distribution of 
infall velocities. In accretion zones with smaller infall ve- 



For TW Hya and V4046 Sgr, where we modelled the accre- 
tion shock in detail, the largest possible post-shock speed 
is 125 km s^^, a quarter of the free-fall velocity. All emis- 
sion redward of this value must originate in the pre-shock 
radiative precursor. 

During the cooling behind the shock front, the velocity 
decreases further. The emission from C ill and O vi peaks at 
very low velocities. For pre-shock velocities of 300 km s'-^, 
O VI is formed immediately behind the accretion shock. Its 
velocity reaches a maximum value of 75 km s~^, and de- 
creases as the cooling proceeds. Depending on the stellar in- 
clination and the spot latitude, the projected velocity may 
be even smaller. The resulting line profile has its maxi- 
mum close to w = km s^^ with a broad wing on the 
red side. Observationally the Ovi hne at 1032 A is best 
suited for analysis, because the other member of the dou- 
blet at 1038 A is more strongly absorbed by H2 and possibly 
by Civ fines. The observed FWHM is about 200 km s^^ 
(Fig. [7]) . Adding turbulent broadening to our 1-D simula- 
tion results in more symmetric line profiles, the red wing 
can be fitted with Wturb = 150 km s~^, but no blue wing is 
observed. In principle, this could be due to wind a bsorption. 
In the case of TW Hya however, iJohns-KruU fc HerczeS 
(I2007D were able to demonstrate that a hot wind is incom- 
patible with, for example, the observed C IV line profiles. 



5.3.3. Model extensions 

Depending on the accretion geometry, radiation from the 
post-shock zone may be reprocessed in the pre-shock ac- 
cretion flow. This will form a photoionised radiative pre- 
cursor, sufficiently hot to emit in Cm and even Ovi lines. 
The emitting region however is located close to the shock, 
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where the material is moving at close to free -fall speed, and 
emission is redshifted by about 500 km s~ MLamzinl ()2003[ ) 
performed detailed 3- level-atom calculations supporting the 
above considerations with numerical simulations. He found 
that the total emission should be observed as two distinct 
peaks, one close to free-fall velocity about 500 km s^^ and 
the other formed post-shock. Only if the line of sight is 
perpendicular to the direction of flow these two peaks can 
coincide, but close to the rest wavelength. For TW Hya, 
which is seen close to pole-on, this would require accretion 
in an equatorial layer. 

Figure [7] shows that in principle a large turbulent broad- 
ening could explain the observed line width, but it requires 
turbulent velocities larger than the post-shock bulk veloc- 
ity. Nevertheless other plasma processes such as magnetic 
waves may supply the required turbulence. 

The line widths could be explained if the emission orig- 
inates in the funnels well before entering the radiative pre- 
cursor, in a region high above the stellar surface. A single 
feature corotating with the star would be observed moving 
through the line profile, as the star rotates, and a combi- 
nation of many emission regions may result in a very broad 
line. In this case an efficient and currently unknown heat- 
ing mechanism in the accretion funnels is needed to form 
the hot ions. The resolution of the current X-ray telescopes 
Chandra and XMM-Newton is insufficient to detect line 
broadening less than about 1000 km s~^ or line shifts less 
than a few 100 km s~^. We are therefore unable to com- 
pare our line-profile measurements with those for hotter ion 
lines, which are most likely accretion-dominated. 

5.4. Turbulence in boundary layers 

Turbulence is an obvious way to generate significant line 
widths. Turbulent flows convert bulk kinetic motion into 
heat and vorticity and a flow along a boundary layer would 
be an obvious origin of turbulence. If, alternatively, ac- 
cretion proceeds in an equatorial boundary layer, the dif- 
ference in velocity is given by the difference between the 
stellar rotation Vrot, and the velocity of material accret- 
ing close to the Keplerian velocity at the stellar surface 

t'Kepier = ^J^^ ~ 440 km s^^, whcrc G is the grav- 
itational constant, and the stellar parameters radius i?* 
and mass M ^, for t he example of TW Hya are taken from 
IWebb et al.l (|1999D . The equatorial velocity can be esti- 
mated from Vea = ^^^^ « 40 km s"^ 

^ sin 2 

In order to estimate the dissipated energy, we use the 
turbulent Crocco-Busemann relation, which is valid for tur- 
bulent flow a long a surface. We follow the treatment by 
IWhitd (|l99ll ) and neglect the initial temperature of the 
infalling flow, since the accretion disk is cool, and the pho- 
tospheric temperature. We obtain a relation between tem- 
perature in the turbulent layer Tturb and velocity u of flow: 

u 

T'turb(M) = -TT—i^V - u) , (3) 
ZCp 

where Av = UKopicr — I'oq, u is the velocity of the gas 
relative to the photosphere and Cp has the usual thermo- 
dynamic meaning as specific heat at constant pressure. 
For an ideal gas, Eqn ^ simplifies to Tturb — 2.3 x 

IO^Kt — — — -T3-(u — u^)i with u = 4^-, so the maximum 

(lOOkms-i) ^ ^" 

temperature is reached where the infalling gas has slowed 



down by ^ . For At; « 400 km s ^ , the maximum temper- 
ature is 9 X 10^ K and 2.3 x 10^ K for half the Keplerian 
velocity. In the example of TW Hya, in both scenarios the 
formation temperature of O vi is reached. Because the line- 
of-sight is perpendicular to equatorial accretion, the peak 
of the velocities should appear close to the rest wavelength 
with a line width comparable to the turbulent velocity. In 
this scenario, it is clearly difficult to explain the shift of the 
line centre, but this process may contribute to the emis- 
sion responsible for a broad component, especially in stars 
with with only slightly shifted emission such as DF Tau, 
V4046 Sgr and TWA 5 or in TW Hya during the GHRS 
and STIS observations. For TWA 5, the scenario further- 
more explains how accretion, as supported by a measured 
Ha EW of 13.4 A, could proceed. This is in spite of the 
fact that the X-ray radiation originates in a low-density re- 
gion, probably a corona, in contrast to that expected in a 
scenario of magnetically-funnelled accretion in hot spots. 

5.5. Red-shifted emission 

All UV observations of TW Hya show consistently red- 
shifted lines. The only possible explanation are matter flows 
onto the stellar surface because the far side of the star is 
blocked from view by the star itself. On the other hand, 
as argued above, the absence of a double-peaked profile 
with pre- and post-shock emission implies that accretion 
spots are unlikely to explain their origin. We therefore now 
discuss a contribution from the a ccretion funnels. Three- 
dimensional MHD simulations of iRomanova et al.l ((lool) 
give accretion spots with inhomogeneous velocity profiles, 
where the flow speed is close to free-fall at the spot centre 
and much slower at the edges. Their simulations use ideal 
MHD equations and do not take into account the heating 
by turbulence caused by the difference in velocity within 
the accretion streams. In general the discussion in Sect. 15.41 
applies with velocity differences of the order of half the 
Keplerian velocity, so significant turbulent heating could 
arise in the funnels themselves, if turbulence is not sup- 
pressed by the strong magnetic fields. Further modelling 
and clarification of the accretion geometry is required to 
check the resulting emission measures. Figure [7] shows that 
the shock, given some turbulent broadening in the post- 
shock zone, can account for a significant proportion of the 
observed emission in accordance with our post-shock cool- 
ing zone simulations. 

Interestingly, in MS stars li nes red-shifted by a 
few km s~^ are also observed ()Avres &: Linskvl Il980l : 
lAvres et al.lll983h . so there has to be an emission mech- 
anism unrelated to accretion, usually attributed to down- 
flows in the stellar transition region. This phenomenon may 
occur in CTTS as wefl. 

5.6. Other emission origins? 

Emission in the disk surrounding the star may in principle 
explain broad M-shaped lines, but th e temperature at th e 
disk surface is only a few thousand K ([Herczeg et al.l[200^ . 
Neither Cm nor Ovi forms there. 

Also, for e.g. TW Hya, where exceptionally broad lines 
are detected, the projected velocity of the disk at the inner 
disk edge is only « 20 km s~'^. We observe no variation in 
the peak shift between the two observations of TW Hya, 
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SO it seems unlikely that it is produced by a single feature 
moving with the rotating disk. 

If the lines originate in an extended chromosphere, from 
the ratio of equatorial velocity and line width we can esti- 
mate the height of the emitting structures and find, using 
the values for Ovi 1032 A from Tabled that the emission 
is placed at the stellar surface for T Tau and at about three 
steUar radii for RU Lup, V4046 Sgr, GM Aur and TW Hya. 
This can only explain line shifts if we assume upflowing ma- 
terial for RU Lup and downflows for TW Hya stable over 
all exposures, which seems unlikely. 

5.7. Variability 

In theoretical simulations accretion, jets and 

wind generation are hig hly ^variable processes 

()von Rekowski fc Brandenburg! |2006| ). Comparing the 
fluxes of TW Hya in the multiplet Cm 1175 A , in 
observations acquired using FUSE and HST, we find that 
the flux decreases by approximately one quarter over three 
years. In contrast, a flux difference of over a factor of five is 
measured for RU Lup observations (Hcrczcg ct al. 2005). 
For TW Hya, there are also FUSE observations separated 
by several years. Unfortunately, the earlier observations 
are very short and of low SNR. Comparing the GHRS and 
STIS observations, separated by 4-8 years, the observed 
line width and shifts appear to be roughly consistent. We 
next considered correlations between these observations 
and FUSE data. Excluding the Cm 1175 A multiplet, for 
which a line-profile analysis is impossible, there are no 
lines that are in common between the two sets of data. It 
is expected however, that lines of Siiv, CiV and Nv will 
originate in emission-line regions for which Cm and Ovi 
emission is detected. In all cases, the line widths fall in the 
range 100-300 km s~\ and for RU Lup, T Tau and DF Tau 
the lines are all consistently blueshifted, but in the HST 
observations to a much lesser extent than in the FUSE 
data, conversely, the TW Hya lines are more redshifted 
in the FUSE data. This hints at a strong variability of 
CTTS, where the dominant UV emission mechanism may 
change between winds and accretion. 

Longer monitoring is required to test if the difference be- 
tween the objects in our sample is due to different phases of 
activity and accretion rate or if there are other fundamental 
differences. This is similar to the situation in X-rays, where 
several luminosity changes are observed, mo st likely due to 
coronal flaring on th e time scale of hours (jKastner et al.l 
l2002t [Giinther et "aL| [2006) . but the evidence is inconclusive 
and some variability might also be due to time- variable ac- 
cretion rates. 



6. Summary 

We have presented the hot ion lines originating from C ill 
and Ovi in all available CTTS spectra observed with 
FUSE. The fitted centroids range from about -170 km s~^ 
to +100 km s~^ and the shifts of the Cm 977 A line and 
the Ovi 1032 A and 1038 A doublet are consistent. Most, 
if not all, lines are superrotationally broadened. The blue- 
shifted lines could originate in a stellar outflow, maybe a 
jet; the red-shifted lines are incompatible with current mod- 
els of magnetospheric accretion. The presence of boundary 



layers, winds and extended chromospheres is inconsistent 
with the observed line profiles. 

Furthermore, the absorbing column densities observed 
in X-rays are incompatible with the optical reddening mea- 
sured assuming a standard gas-to-dust ratios and interstel- 
lar extinction law. This situation holds for the most heavily- 
vefled objects AA Tau, DF Tau, T Tau, SU Aur and, al- 
though for lower absolute values of reddening, RU Lup. In 
cases where data are available from both UV and X-ray, the 
corresponding objects demonstrate blue-shifted emission in 
the FUSE observations. We interpret this as a sign of a 
dust-free absorber, which may consist of shock-heated jet 
material. In the ratio between the O vii luminosity and the 
O VI luminosity, all CTTS display a clear excess of soft X- 
ray emission in the O vii lines. Together with X-ray results, 
this illustrates that the signatures of a hot accretion spot 
are most evident in the temperature range 1-2 MK. 
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